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ABSTRACT 


Previous work in the field of water supply contamination from fallout due to 
nuclear attack has indicated that the level of activity brought into the 
water supply by surface runoff from the watershed increases the activity con¬ 
siderably. These studies employed a constant value for the "runoff coefficient" 
to obtain estimates of the degree of contamination resulting from the soluble 
fraction of fallout activity deposited on the watershed. The present report 
is part of an investigation to better define the relationship between rainfall 
and the amount of activity that reaches the water supply due to runoff. It 
concentrates mainly on two areas: the hydrology of surface runoff and ion 
exchange phenomena. 

A detailed qualitative discussion of the relationships between the various 
hydrological parameters is presented. Some of the difficulties encountered 
in trying to predict runoff as a fixed percentage of rainfall are discussed. 

Some data on ion exchange properties of soils with respect to the biologically 
important radionuclides are given. In addition the report induces an 
indication of how the information presented in the report may eventually be 
included in a consideration of the general problem of transport of activity, a 
discussion of uptake of activity by biota and some preliminary calculations 
of water supply contamination for Providence, R.I., one of the cities in the 
OCD "5-City Study". 


The maximum levels of selected contaminants at the Providence water intake from 
a total 30 MT weapons attack were found to be approximately as follows 
(in pc per mi!): 



3r-d9 

Sr-90 

Ru-106 

1-131 

Cs-137 

Ba-l40 

Surface cotroaminat. i only 

5 x 10“ 3 

3 x 10 

4x10 '■' 4 

4x10 ~ 2 

2 x 10" 5 

-2 

3 x 10 

Contamination,incl.runoff 

7 x 10“ 2 

7 x 10 _1+ 

1 x 10 

1 

4X10" 11 

6 x 10“ 1 


These figures indicate ■‘•hat the water supply contamination problem may be 
increased by a factor of 10 to 25 when surface runoff from the watershed is 
considered. 
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I. INTRODUCTION 


Hie general purpose of this study is to evaluate and suninarize the problem of 
water contamination by radioactive fallout in the event of nuclear war. The 
level of fallout that might result from a possible nuclear attack was based 
on current theories of the formation and distribution of fallout. 

The potential degree of water contamination from fallout was evaluated and 
estimates of the biological hazards associated with the ingestion of fallout 
contaminated water following hypothetical nuclear attack presented (l). The 
potability and biological effect of contaminated water supplies were questioned. 

The availability and use of water supply are most important to a conmunity 
during the postattack period. From previous investigations it appears that 
surface runoff may contribute significantly to the contamination of water supplies, 
depending on environmental factors including watershed characteristics 
and meteorological aspects (*). Therefore, it was proposed to conduct a 
detailed analysis of the redistribution of fallout and the contribution 
to contamination from surface runoff during the early postattack period from 
land areas to streams and from watersheds into reservoirs. The analysis 
of the redistribution phenomenon is equally important in and around urban 
centers wh-re the activity concentration may change significantly following 
precipitation, especially for the soluble nuclides. From the previous studies 
it appears that surface runoff from rainfall during the first 2b to bd hours 
following nuclear attack will be critical as far as water contamination is 
concerned. 


(*)fesults of these studies indicated that the contribution to water con¬ 
tamination due to surface runoff may be significantly greater than that 
due to direct contamination by a factor of 10 to 25. When the effect of 
watershed runoff was included, contamination levels exceeding 10”^ (i c/m£ 
for 1-131 and Ba-l^K) were obtained. 
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II. OBJECTIVES 


The objective of these investigations are: (l) to conduct a detailed analysis 
of the degree of contamination of water supply systems due to runoff, (2) to 
use these results to estimate the dose absorbed by human body over different 
periods of water consumption, and (3) to contribute to the OCD 5-city study 
by evaluating fallout contamination of selected water supplies. 

To determine the mechanism by which fallout activity is transported by the 
runoff and contributed to water supply contamination it is necessary to 
study e variables affecting runoff; including hydrologic, meteorologie 
and gf ^ogic aspects, investigate the solubility, ion-exchange, biota 
uptake, transport and redistribution of radionuclides in water and over land 
and to examine the effect of different watershed vegetation and foliage. 

To estimate the potential hazard to the population from the ingestion of 
fallout contaminated water it is necessary to determine the atom concentra¬ 
tion in drinking water. The resultant body burden is calculated for different 
ingestion starting times and over different periods of consumption based on 
a consideration of physical, chemical, biological and physiological factors 
in their relation to water consumption and their interrelationship within 
the human body. 

Io perform computations on water contamination for selected water supplies 
as required by the sponsor to contribute to the OCD 5-city study, the 
previous methods of evaluation are being refined. Existing computer 
programs need to be streamlined to save valuable computer time and new 
computer programs for total atom concentration at any location from a single 
weapon and for multiple weapon attacks are being developed. 
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III. SCOPE 


From the previous work it appeared that the degree of radioactive contamination 
from surface runoff, depending on environmental factors and especially watershed 
characteristics including hydrologic and physico-chemical aspects, may he a 
major factor in the rehabilitation and use of public water supplies. 

Since empirical hydrologic expressions can not be applied directly to problems 
involving both hydrologic aspects and transport of activity, it was found that 
it is necessary to establisn relationships between hydrologic and physico¬ 
chemical pehnomena, such as ion exchange, diffusion, particle size, and other 
phenomena relevant to transport (2)(3)(^)(5) and to formulate these for quanti¬ 
tative use. 

The fallout contamination of water supplies due to surface runoff during the 
postattack period may be evaluated by consideration of the following objectives: 

A. Hydrology and Surface Runoff 

This phase of the project essentially consists of a comprehensive delinea¬ 
tion of the hydrologic cycle leading into methods of predicting runoff from 
rainfall and from antecedent conditions. It is necessary to closely examine 
the variables affecting runoff, such as: (a) slope, (b) ground roughness, 

(c) soil moisture, (d) antecedent rain, (e) terrain, (f) vegetation, (g) season, 
(h) soil composition and/or rock formation(s), (i) wind velocity and direction, 
(j) rainfall intensity, (k) rainfall frequency and duration and (l) watershed 
management and cultivation. 

An analysis of these variables is necessary to obtain quantitative results 
of runoff for the various conditions that will be encountered postattack. 

B. Chemical and Ion Exchange Properties of Fallout Nuclides 

The solubility of radionuclides from fallout in surface runoff will be a 
major path of water contamination in the period after a nuclear war. When new 
ions from fallout particles are introduced into the soil, they compete with and 
replace other ions on exchange sites in the soil. In some reactions with the 
soil the new ions become nonexchangeabie and thus some portion of the radio¬ 
element is rendered unavailable for uptake. Therefore, the types of interactions 
that occur between the soluble radionuclide and the s--.ll constituents determine 
the availability of the radionuclide for uptake from the soil. Similar consider¬ 
ations enter the availability of a given radioelement for uptake by stream 
sediments and biota. 

It is necessary to study the extent of surface penetration of radionuclides 
by natural processes and the extent to which cultivation practices determine 
the initial concentration in the soil, i.e. soil uptake contamination factor 
relationships (5). In addition, it is necessary to investigate the following: 

(a) factors affecting the removal of various isotopes from different types 
of surfaces, including various soil cavers, foliage, vegetation and crops, 

(b) redistribution of various radionuclides alone or in chemical combinations 
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which may be expected following a nuclear attack and (c) ion-exchange properties 
and capacities of various soils, vegetation, stream sediment and biota to 
predict the retention behavior of radionuclides. 

Although results are not always comparable, available data indicate that 
the phenomena of ion-exchange, chemical properties, uptake by biota and particle 
size exert considerable influence cn the amount of nuclides ultimately available 
for surface transport by runoff. 

C. Transport of Nuclide Activity 

To calculate the effect of surface runoff on contamination it is necessary 
to know the soluble fraction of the activity of each biologically important 
nuclide. A study on the purely theoretical aspects of the problem of activity 
transport, in conjunction with ion-exchange, retention by solids, uptake by 
biota and redistribution phenomena has been initiated. As a mathematical 
transport model can be established, problems of uptake, ion-exchange and 
redistribution may then be superimposed to yield the ultimately desired 
results. 

To analyze the transport problem it is necessary to examine available 
information in (a) activity transport by streams, including tracer studies (6) 
and sediment transport (2)(U), (b) overland transport, including surface 
transport (7) and (c) mixing of activity in reservoirs (6). 

The analysis of the redistribution problem is expected to result in a 
mathematical model employing partial differential equations. The transport 
of activity by stream flow involves the theory of turbulent diffusion. As 
water from watershed runoff and streamflow approaches a reservoir or other 
area of surface storage, the initially deposited radioactivity undergoes 
diffusion. The boundary value problems of combined subsurface & overland 
flow, uptake of activity by soil and transport of activity by streams should 
be programmed for digital computer solution. 

D. Estimation of Absorbed Dose 


Following the evaluation of the concentration of water contamination 
from surface runoff, the potential short and long term hazards and radio¬ 
biological effects will be estimated. In connection with this analysis, the 
specific methods of calculating contamination levels will be reviewed and 
wherever possible the previous methods of evaluation refined. 

To evaluate water contamination due to surface blasts of 5, 10 and 20 MT an 
arduous procedure had been followed. The difficulty was due to the lack of 
a successful computer program to evaluate nuclide solubility contour ratios. 

By incorporating the recently made available contour ratio values (N/l) into 
the existing progra® Tor intensity evaluation (i), the number of soluble atoms 
per sq ft (N) are found directly. This interpolation subroutine for computing 
the atoms/sq ft, N, directly will involve modifications of the program entitled 
"Estimated Intensity at Any Location in the Fallout Region for a Single Weapon 
Yield" (8). The program for multiple weapon yields will be modified similarly; 
it will give a read-out in atoms per unit area. The multiple weapon program 
previously presented (l), while efficient for two weapons, has proven to be 
quite uneconomical and too complicated for the poly-weapon situation. In 



addition only the total intensity at a point was calculated. The new multiple 
weapon program is expected to give tne total atom concentration for the entire 
watershed as well as for any point within the area of the attack. This latter 
output will be of special value to the evaluation of contamination due to runoff. 
For these reasons, a more sophisticated program is required for our specific 
needs and is currently under development. The new program will employ various 
subroutines and reserve the MAIN program for direct computation and summation 
of atoms per unit area. 

With the aid of these computer programs it will be possible to estimate 
the absorbed dose in the human body over different periods of water consumption 
from a single surface burst and to maximize the contamination level due to 
several bursts and saturation raids. 



’f 
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IV. HYDROLOGY AND SURFACE RUNOFF 


A. Introduction 

The science of hydrology treats the behavior of water in the atmosphere, 
on the ground and underground. Although always important to the survival 
of man, it is only in comparatively recent years that attempts have been 
made to improve the empirical relationships which have governed this field. 

One of the primary aims of the science of hydrology has been the study of 
the relationship between rainfall and runoff. 

In the previous report (l) the contribution of fallout on the watershed 
to water supply contamination was estimated by the use of an approximate 
runoff coefficient in the so-called Rational Formula, Q = CIA, for runoff. 

This estimate was used because it was not possible to investigate the problem 
to any greater depth in the time available. However, the estimate indicated 
that this contribution to contamination, namely, the amount brought into the 
water supply by surface runoff following a rainstorm, may in seme instances 
be the largest contribution to the contamination. 

The central purpose of this project is to improve the reliability of the 
estimate of the water supply contamination brought into the water supply by 
surface runoff. The physical phenomena that are relevant to surface runoff 
are so complex that little far reaching theoretical work has appeared in the 
literature. In fact, even the transport equation for water in a defined 
channel has been solved analytically for only a few special cases (2). 

The field of hydrology has been confined for the most part to empirical 
relations with "coefficients" and "capacities", some of which are not well 
defined Jusofar that the definition leads to a means of measuring the quantity. 
Furthermore, the relevant parameters are not always explained with the 
presentation of data. 

The hydrologic cycle, the movement of water as it relates to the earth, 
is an extremely complex system. To solve the problem of contamination from 
surface runoff requires a knowledge of (a) the fraction of activity due to 
specific radionuclides which are soluble in the runoff, (b) the speed with which 
the soluble portion of the activity may be transported overland by runoff to 
streams, lakes and reservoirs, and (c) the exchange reactions and other 
physico-chemical phenomena that affect the amount of activity that will 
actually reach the water supply. A series of modified transport equations 
may be ultimately used to solve this problem analytically. However, before 
any part of the overall problem can be solved it will be necessary to under¬ 
stand the dynamic behp lor of the hydrologic cycle, in particular the runoff 
phase, in order that one prediction of the amount of runoff from a particular 
storm may be possible. 

B. Fixed Limits of the Hydrologic Cycle 

The oceans contain 97$ of all the water in the world, or one quadrillion 
(1015) acre-feet. The total amount of fresh water on the other hand is 


''•»*«*- • — ■- ^ tl . - ^itrftffii-rttrirv aiif - ,r-' nif -urm-nit -- ■ -iT‘- .mss «<»-*• 
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32 

estimated at 33 trillion (33 x 10 ) acre-feet. It is distributed roughly as 
shown in Table I. 

TABLE I 


Worldwide Distribution of Fresh Water 


Source 


Quantity ($) 


Polar ice and glaciers 75 

Ground water 

At depths < 2,500 ft 11 

At depths > 2,500 ft, 

but < 12 ,500 ft- 14 


Lakes 0.3 

Rivers 0.03 

Soil (moisture) 0.06 

Atmosphere 0.035 


These are estimates of a static relation, and whereas the water content 
of the atmosphere is relatively small at any given time, immense quantities 
of water pass through it annually (9). For example, the annual precipitation 
on land surfaces alone is 7.7 times as great as the moisture contained in the 
entire atmosphere at any one timej that is, about 30 times as great as the 
moisture in the air over the land (10). 

As can be seen from the figures above, an understanding of the "flow" 
nature of the cycle is necessary to appreciate the tremendous quantities of 
water in transit through the atmospheric and continental phases of the cycle. 


C. The Hydrologic Cycle 

Due to the complexity of the hydrologic cycle, water particles may be 
considered as following one of three paths. Of the water that reaches the 
land surface by precipitation, some may evaporate where it falls, some may 
infiltrate into the soil, and some may run-off overland to evaporate or 
infiltrate elsewhere or to enter streams (ll). The larger part of the 
water falling on the land surface passes back to the air, either directly 
by evaporation or through the bodies of animals and plants in respiration 
and in transpiration. This evapo-transpiration accounts for 71$ of the 
five billion acre-feet of rainfall that fall over the continental United 
States eauh year. 

A smaller portion, the remaining 29$ of the water reaching the land 
surface runs off and either sinks into the ground or stays on the surface 
until gathered by rivers that carry it back to the sea. Of the water that 
infiltrates into the ground, some may be evaporated, some may be absorbed 
by plant roots and then transpired, and some may percolate downward to ground 
water reservoirs (12). Water from these underground reservoirs may move 
laterally until it is close enough to the surface to be subject to evapora¬ 
tion or transpiration, some may reach the land surface and form springs, 
seeps, or lakes, and some may flow directly inuo streams or into the oceans. 



•» .r •- * «'•£»• V 
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Of the water that remains on the surface, some may accumulate in lakes 
and surface reservoirs, sane may be lost through evaporation or transpira¬ 
tion or riparian vegetation, some may seep downward into ground water reservoirs, 
and some may continue on to a saline lake or to the ocean (ll). A very small 
fraction is taken up in the bodies of plants and animals. Plants use its 
contained hydrogen in forming their tissues. The hydrogen is later recombined 
with oxygen through animal and plant respiration; the water thus produced 
returning to the air. The hydrologic cycle is completed by evaporation from 
the earth's saline water bodies and circulation of water vapor in the atmosphere. 

D. Hydrologic Disposition of Precipitation 

Rain and snowmelt seep into underground water-bearing strata (aquifers) 
and run off saturated-soil surfaces into streams, lakes, and rivers. Of 
this total precipitation, 70 to 75$ is returned to the atmosphere by 
evapotranspiration from plants and vegetation, and by direct evaporation 
from moist soils, lakes and streams. This water is not available for man's 
volitional use. Slightly more than half of it is beneficial, however, as 
it sustains the growth of forests, grasslands and non-irrigated farmlands. 

The portion of precipitation that seeps into underground aquifers is 
particularly important, because the gradual overflow from full aquifers 
accounts for most of the regular streamflow (as distinct from that resulting 
from surface runoff). 

1. Precipitation 

Precipitation results when water vapor in the atmosphere is cooled 
sufficiently for some of the moisture to condense on dust, particles or other 
condensation nuclides to form water droplets. Natural cooling of air 
masses results from (a) solar radiation or (b) lifting. Cooling by radiation - 
the cause of dew and frost - is relatively unimportant. However, cooling by 
lifting is the primary cause of precipitation, and may result (a) from topo¬ 
graphic conditions, (b) from the convergence of weather fronts, or (c) from 
thermal convection. 

In this discussion precipitation was assumed to be rain falling 
on unfrozen ground. The presence of snow, ice, or fro^f- may have a great 
effect on the disposition of water in the hydrologic cycle. Interception 
losses may be much greater if the precipitation is snow. Evaporation losses 
from snow surfaces may also be very rapid. Much water may be stored temporarily 
on the ground until released by melt or evaporation. Frost in the ground may 
either increase or decrease the infiltration capacity, depending on the 
moisture content of the soil; a frozen soil with nigh moisture content may 
approach impermeability. A thick blanket of snow will prevent freezing in the 
ground, and snow-melt at the soil surface may keep the soil near saturation. 

Freshly fallen snow is relatively "dry", i.e., of low density. The 
absorption of heat, even at temperatures below freezing, increases the density. 
Snow-melt is insignificant until maximum density is reached, then additional 
heat causes relatively rapid melting. Considerable rainfall can be absorbed 
by dry snow, but even a small amount of rain when the snow is close to maximum 
density, may result in complete sncw-melt and in runoff exceeding the precipitation. 
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The rate of melting is proportional to the heat absorbed and, for this reason, 
is a gradual process unless'hastened by rain. Although the ultimate disposition 
of sncw-melt is similar to that for rain, the initial losses of rainstorms, 
consisting of interception, wetting of surface, and depression storage, are 
not so important in snow-melting. 

2. Interception 

Precipitation may be retained in foliage and evaporated or it may 
temporarily collect on leaves or needles, then drip or run down stems or 
crowns to the ground. That which is caught in the foliage and evaporated 
contributed nothing to water yields. The amount intercepted depends largely 
upon intensity, amount, and kind of precipitation and the density, type and 
height of vegetation. 

Interception losses are greatest for areas of coniferous forest in 
winter, but rain interception may be appreciable for any type of vegetation, 
including grasses. The degree of interception differs according to whether 
the vegetation is coniferous, broadleaved, evergreen, or deciduous, has large 
or small crowns, short or long branches and rough or smooth bark. A light 
showerfmay be almost entirely intercepted. As the amount which can be held 
by leaf surfaces is limited, interception becomes less important with prolonged 
or heavy rains. Thus, only a small proportion of a heavy rain is intercepted; 
while snow is more readily held than any other form of precipitation. 

Evaporation from leaf and other surfaces of vegetation probably is 
slight during a storm, because of a high relative humidity, but most intercep¬ 
tion loss is evaporated within a short period after precipitation ends. 
Practically no moisture is absorbed by plants from wet leaf surfaces; the 
immediate cooling effect of a summer shower on vegetation is generally caused 
ly a reduction in transpiration rate. 

Interception losses range from nearly 100$ for light showers to about 
25$ for heavy rains. The total annual interception over the United States 
has been estimated as high as M}$ of the annual precipitation in mature stands 
of coniferous forests and 15$ of the annual precipitation as a general 
average for all forests (13)• However, conditions vary so greatly even in 
small areas that general quantitative statements are of questionable value. 

3. Infiltration and Percolation 


Precipitation may (a) fall on a water surface, (b) be intercepted by 
vegetation, or (c) fall directly on the soil. A particle of water reaching the 
soil may (a) be returned directly to the sky by evaporation, (b) flow overland 
toward a stream, or (c) enter the ground b;y infiltration. The process of water 
entering the soil is called infiltration, and the movement of water to deeper 
depths in the soil is called percolation. More specifically, percolation is 
the lovement of water under hydrostatic pressure through a rock or soil, 
excluding turbulent flow through large openings (l4). 

Evaporation from land or water surfaces during a storm is small 
and surface runoff does not begin until the soil surface is thoroughly 
wet and shallow depressions are filled; infiltration, however, begins 
with the first drops of rain. Water entering the soil first replenishes 
moisture deficiencies in the soil, and then contributes to streamflcw and 
ground-water supplies. 
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Infiltration is determined by conditions at the immediate surface, 
whereas percolation is controlled by subsurface conditions; the rate of 
both often begins relatively high and decreases to a much lower and more 
or less constant rate as the rain continues. Molecular attraction is 
effective in pulling water into the soil until the upper part of the soil 
becomes wet. However, after this occurs, gravity is the principal force 
at work, and infiltration becomes a function of the permeability of the 
soil. Consequently, soils having coarse texture or of large aggregates, a 
layer of litter and humus on the surface, insects and small animals working 
in the soil and abundant roots of vegetation to make spaces for water move¬ 
ment and storage, have high infiltration and percolation rates. Other 
conditions that may cause variations of infiltration capacity are: (a) soil 
moisture content; (b) state of cultivation; (c) perforations of the surface 
soil and subsoil, such as those left by earthworms and decayed roots; (d) 
packing of the soil surface and the clogging of the soil pores by i.ne particles 
washed down from the surface by rain; (e) temperature changes; (f) shrinking 
and swelling of surface soils, including sun-checking during dry periods; 
and (g) depth to less permeable strata (14). 

Since many of the above factors cannot be duplicated in the laboratory, 
it is usually necessary to measure the infiltration capacity of soils in 
place. Among the methods are the following: (a) measurement of the rate of 
intake of water on very small areas bounded by metal rings or tubes; 

(b) measurement cf the rate of intake of water in areas artifically flooded 
by irrigation; (c) measurement of runoff of water applied to small areas by 
rainfall simulators; (d) comparisons of measured precipitation with measured 
surface runoff; (e) lysimeter studies; and (f) detailed measurements of soil 
moisture content at various depths. 

Soil, or field moisture, is the moisture held by the soil close to 
the surface of the ground, usually within about three feet. Field moisture 
is held by capillarity, i.e. trapped between soil grains or held in colloidal 
gels far above the zone of capillary water that overlies the water table. 

This moisture supports most plant life and permits biological activities in 
the topsoil to continue by acting as a storage reservoir, replenished during 
rain and drawn upon during dry spells. 

The most important factor that affects infiltration is the condition 
of the soil pores. Soil structure, especially near the surface, is extremely 
variable because so many different things act to change it. The most stable 
soils are sands, which support little plant life, and which undergo little 
change with the addition of water. Less stable are clays whose colloidal 
constituents swell with the addition of water. 

The range of infiltration capacities for various soils has been found 
to be quite large, and quantitative measurements are applicable only to the 
soils and conditions at the time of measurement. Initial infiltration 
capacities may exceed 10 in./hr. and the ultimate or nearly constant capacity 
reached toward the end of a rain may approach or reach zero. Initial capacities 
are difficult to estimate, as the capacity usually drops rapidly after rain 
starts. The ultimate capacities may be measured more accurately but may vary 
considerably, even for the same soil,. 
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Horton ( 15 ) introduced the concept of infiltration, defining infiltra¬ 
tion capacity, f, as the maximum rate at which a soil in any given condition 
is capable of absorbing water. One of the most striking characteristics of 
infiltration capacity is its extreme variability with respect to both area 
and time. At the beginning of any storm the infiltration capacity, f 0 , is 
likely to be high; it thi i decreases rapidly during the first hour or so and 
finally levels off and approaches a constant value, f c , after the next hour or 
two. The ratio of f 0 to f c depends greatly upon permeability. For a well- 
compacted clay it is high, and for a coarse, sandy soil it is relatively lew. 

This characteristic may be more readily understood when it is realized that 
the infiltration capacity of any particular area is determined at the ground 
surface, ordinarily it is the nature of the openings that exist in the top 
£ in. of the soil that determines infiltration capacity. However, where 
subsurface storm flow occurs it is the condition of the surface layers of 
the relatively impervious strata which determines the infiltration capacity. 
Inasmuch as in all cases f is dependent upon the character and condition 
of this layer, it is readily ' "derstandable that any disturbance of that layer 
may completely change the in ation capacity. From all this it appears 

that infiltration capacity is not a permanent characteristic of a watershed 
that is comparable with size, length, and similar fixed properties. Instead, 
where the objective is the determination of surface runoff, it is necessary 
to determine the range within which the average infiltration capacity of a 
basin varies from time to time; then knowing the factors that cause those 
variations and knowing the influence of each, estimate the value that should 
prevail at any given time. Thus, with the miniroam, maximum, or average 
infiltration capacity existing during any design storm, it is possible to 
determine the greatest, smallest, or average runoff that will result from 
such a storm. 

The typical infiltration-capacity curve starts with a high value, 
i'o, because of the initial soil conditions, then dropping rapidly during 
the early stages of the storm and finally levelling off and approaching a 
constant value, f c . For different watersheds this ratio, f 0 /fc> has a 
wide range, depending on the amount of interception and depression storage, 
the type, texture and condition of the soil as previously mentioned. 

Horton (l6) suggested that it would decrease exponentially with time during a 
storm and approach a constant rate, usually after a period of 1 to 3 hr. and 
that it may be represented by an equation of the form: 

f = f c + (f o -f 0 )e" kt (1) 

where: f = infiltration capacity at any time (inches/hour) 

f c = constant, or ultimate, infiltration capacity 
(inches/hour) 

f 0 = initial infiltration capacity (inches/hour) 
e = Jfapierian' base 

k = a decay constant for a given curve 
t = time from the beginning of precipitation (hours) 
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Attention should he called to the fact that the infiltration-capacity 
curve is, in practice, not an f curve but, rather, it is an average infiltra¬ 
tion capacity, f a , curve. In other words it dees not represent instantaneous 
values of infiltration capacity as it existed throughout the storm period, 
but instead it represents the average infiltration capacity for each of the 
several periods of high storm intensity. The infiltration approach to runoff 
is based on the use of infiltration indices. The B-index. for example, is 
the average rainfall intensity above which the total mass of rain equals the 
total mass of the observed runoff. The remainder of the total precipitation 
consists of gross basin recharge, i.e. surface retention plus infiltration. 
Another concept, the W-index, has been defined as the average rate of 
infiltration, usually measured in inches depth per hour for a given area, 
during the time rainfall intensity exceeds infiltration capacity: 

F 

W = -fjr = j (P - Q.g - s fi ) (2) 


where: W = average rate of infiltration, when rainfall intensity 

exceeds infiltration capacity (inches/hour) 

= total amount of infiltration (inches) 

T = time during which rainfall intensity exceeds infiltration 
capacity (hours) 


P = precipitation (inches) 


Q s = observed surface runoff from the storm (inches) 
S e = total surface retention (inches) 


This approach has been successful in estimating maximum flood flows 
which consist almost entirely of surface runoff, when the initial moisture 
condition of the soil is quite uniform. The W-index is estimated and an 
infiltration curve superimposed on the rainfall plot. The surface retention 
is then estimated and added to the infiltration curve. The surface runoff 
is represented by the area between the rainfall curve and the infiltration 
plus surface retention curve (17)- 


4. Surface Runoff and Subsurface Storm Flow 


As long as the rate at which rainfall reaches the soil surface is 
less than the infiltration capacity, all the available supply of water sinks 
into the soil. As rain continues, plant surfaces become saturated, the 
interception-loss rate declines, and the infiltration capacity also decreases. 
Whenever the supply rate of rain exceeds the infiltration capacity, sha-low 
depressions begin to fill with water, and when these depressions are filled 
to overflow level, water begins to move by overland flow toward streams. 
Overland flow usually reaches a tiny rivulet or channel within a short distance. 
The time of travel from raindrop to stream channel depends on distance, slopes, 
and surface conditions end may also be affected by depth of flow. Generally 
the time of travel is a matter of minutes. 
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"Surface runoff" is the water which reaches a stircm by overland flow 
and is a residual equal to precipitation minus the total evapo-transpiration 
losses and the total ground-water flow. It can be directly measured cm small 
sample plots and cm ephemeral streams, but cn larger streams the hydrograph 
of streamflcw is complicated by ground-water inflow and channel storage. In 
this area, the analysis and synthesis of the hydrograph is one of the funda¬ 
mental problems of stream hydrology. 


A division does not take place at the soil surface between soil 
moisture and surface runoff. As indicated above, not all overland flew 
reaches a stream channel - thus, the distinction between overland flow and 
surface runoff. This water which flows through the soil but above the 
water table, is called "subsurface storm flcv." Under seme conditions, 
when a relatively impermeable layer retards or prevents the percolation of 
water downward and diverts it back to the surface or into stream channels, 
subsurface storm flow may reach stream channels almost as rapidly as doe3 overland 
flow. However, if the soil above the water table becomes saturated, it 
constitutes a temporary or perched water table and may contribute a substantial 
amount of water to streamflcw during or after a storm. A sharp division between 
surface runoff and ground-water seepage flow is therefore, impossible. 

5. Streamflcw 

Streamflcw is measured by recording the stage or elevation of the 
water surface at a given station. The actual average flow velocity in 
the cross-section at the station is measured with a current meter for a 
number cf flow conditions at different stages. These values are used to 
get a stage-discharge relationship, or rating curve. The discharge Q is 
computed as: 

Q z A i v i 

where; A.^ = cross-sectional area of the part of a stream at a gaging 
station, and 

= average velocity in that part of the stream. 

The average velocity, V, of the stream at a stauion therefore is: 

£ A v 

TT *“■ “ . 


The discharge may be determined by the slope of the energy line 
(approximated by the slope of the water surface) as follows: 


Q, = A R 2 /-’ s 1 / 2 


where: Q = discharge (cfs) 


A = cross-sectional area of flow (sq ft) 

R = hydraulic radius: ratio of A to contact line of A witi. 
river bed and banks 

S = slope 

n = Mannings roughness coefficient 
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For two discharges at the same station and stage: 




(4) 


6. Transpiration 

The biological sciences of botany, ecology, and plant physiology are all 
concerned with this phase of the hydrologic y cycle. Plants absorb water, 
principally through their roots, use the water in several involved ways in 
their physiologic processes, and give off water to the atmosphere largely 
through water-vapor diffusion through the pores or stomata in the process 
called "transpiration". The amount of water held in storage by a plant is 
less than 1# of that lost by it during the growing season (l4). Therefore, 
from the hydrologic standpoint plants are pumps which remove water from the 
ground and raise it to the atmosphere. 

In some areas where seasonal drought is usual, the amount of water 
removed from the soil by transpiring plants depends largely upon the depth 
of root penetration. However, since transpiration is essentially a process 
of evaporation, it depends largely on the same factors that influence 
evaporation from land and water surfaces - namely, solar radiation, temperature, 
wind, and relative humidity. 

.Although numerous experimental methods of determining transpiration 
rates have been developed, their general applicability to hydrologic studies 
is varied, since many of them were developed for the study of specific 
problems. One method consists of weighing freshly cut parts of plants, 
immediately after cutting and periodically thereafter until wilting starts. 

It is based on the assumption that transpiration continues at the normal 
rate immediately after cutting. Another method is by potometer measurements. 

A pocometer is a vessel containing water into which the cut end of a plant 
or leaf is inserted. After sealing, measurements are made of the amount of 
water removed from the vessel. A third method is by phytometers. These differ 
from potometers, in that they contain soil in which the whole plant is grown, 
thus approaching natural conditions. The closed phytometer is used extensively. 

The quantity of water involved in the transpiration process is 
important as can be seen from the following: "An acre of corn gives off 
about 3>000 to 4,000 gallons of water each day while a large oak tree gives 
off about 40,000 gallons per year" (l3). 

7. Evaporation 

Evaporation from soil surfaces varies roughly in the same manner as 
does transpiration and usually can not be separated from transpiration 
losses. Soil evaporation rates vary within wide limits, from approximately 
the maximum rate for free-water surfaces to zero (l4). Vegetation shades the 
soil and reduces the soil evaporation, but transpiration usually exceeds this 
reduction, so that plants increase the total losses. 
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Because of tht nearly insurmountable difficulties of making separate 
measurements of the two phenomena, evaporation and transpiration are generally 
lumped as "evapo-transpiration." If we neglect storage in the soil moisture 
and in ground water, then the difference between total precipitation and total 
stream flow is the total evapo-transpiration. 

The measurement of evaporation from free water surfaces by "direct means" 
requires the measurement of all inflow, outflow, and scorage, with the difference 
assumed to be evaporation. There are few natural water surfaces for which 
evaporation losses can be adequately measured by this method. 

There are several approaches to indirect measurement of evaporation, 
but the various factors affecting it are difficult to estimate. Some of 
these factors are temperature of the air and water, differences in vapor 
pressure, humidity of the air, solar radiation, wind movement, barometric 
pressure, and chemical quality the water. 

The maximum possible evaporation rate has been given by Meyer (19): 


dE 

dt 


< e s - e a> (1 * 


-) 

10 ' 


(5) 


where: 

— = maximum evaporation rate (inches/day) 

e = saturation vapor pressure (mb or in. of Hg) 
s 

e = existing vapor pressure in air (mb or in. of Kg) 

8 

W - wind velocity about 25 ft. above the surface (mph) 

C - a numerical coefficient ranging from O.36 for ordinary lakes 
of about 25 ft. depth to 0.50 for wet soil surfaces, small 
puddles and shallow evaporation pans. 

Barometric pressure has only a small effect on evaporation, A 
drop of the barometric pressure from 30 in. to 20 in. increases the 
evaporation rate by 20$. 

Water pollution decreases the rate in proportion to the percentage 
of dissolved solids, i.e., sea water with 35>000 ppm of total solids 
(96.5$ water) evaporates 9^.5$ as rapidly as fresh water. 

When considering evaporation from land, the availability of the 
water is a factor. Evaporation opportunity is 100$ for lakes and streams. 

a. Reservoir Evaporation 

The principle of conservation of mass as expressed in the water 
budget equation may be used to determine evaporation from a reservoir: 

p + i = e + ~£- + o 

A 


(6) 
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precipitation on reservoir 

inflow 

outflow 

evaporation 

change in storage per unit area of the reservoir. 


This equation applies either instantaneously in which case the 
qualities are to be considered rates or over a given period of time in 
which case the quantities are to be considered accumulated totals. There 
may be a large error in the measurement of evaporation by this method, 
especially over a short period of time, due to the inherent difficulty 
in measuring the subsurface components of I and 0. 

The most common method is by means of evaporation pans after 
establishing a relationship between pan evaporation and reservoir evaporation. 

If the net evaporation is greater than precipitation, streamflcw 
always decreases as a result of the construction of a reservoir, and 
usually, but not necessarily vice-versa (17). 

b. Monthly Evaporation 

A graphical relationship may be developed, similar to Meyer's 
equation for daily evaporation, to give monthly evaporation (20), as follows: 

E = (7.8 + 0.32 v w ) (e g - e a ) (7) 

where: E = monthly evaporation (inches) 

v^ = wind velocity (mph) 

e = saturation vapor pressure (inches of Hg) 
s 

e & = actual vapor pressure (inches of Hg) 

E. Hydrographs 

The hydrograph is a plot of discharge or runoff as a function of time. 

The study of the hydrograph of a river before, during and after a storm 
is useful in determining how much of the storm precipitation reaches the 
river as direct surface runoff or storm runoff. The study of the hydrographs 
of a series of storms of various durations and under different conditions 
is a first step in establishing a rainfall-runoff relation for a given basin. 

1. Hydrograph Composition 

The four components of runoff are surface runoff, interflow, ground- 
water flow and channel precipitation. The most relevant parameters and 
their influence on these four components of runoff are: rainfall intensity 


where: P = 

I » 

0 = 

E = 

AS _ 
A av~ 
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and duration, infiltration rate and volume, and soil moisture deficiency. 
Generally speaking, four different shapes of hydrographs occur depending on 
the relative magnitude of these parameters: 

Case I : Rainfall intensity, i < infiltration rate, f. ; and volume 
of infiltrated water, F i < soil moisture deficiency 

Because of the second assumption, there will be no contribution 
to runoff from either interflow or ground water flow from this storm. 

Because of the first assumption, there will be no surface runoff. Therefore 
the only runoff contribution will be that due to rainfall directly on the 
stream (channel precipitation). The hydrograph for this case will show 
only a slight increase of the discharge with respect to time over the 
expected base flow. 

Case II : Rainfall intensity, i < infiltration rate, f , and 
infiltration volume, > soil moisture deficiency 

In this case, after soil moisture has reached its maximum value, 
there is an additional contribution of interflow and ground-water flow 
to streamflow. 

Case III : Rainfall intensity, i > infiltration rate, f., and 
infiltration volume, F^ < soil moisture deficiency 

In this case there are contributions to streamflow from surface runoff 
and from channel precipitation, but, due to the second assumption, no 
additional contribution to the existing base flow from interflow or ground- 
water flow. 

Case IV ; Rainfall intensity, i > infiltration rate, f., and 
infiltration volume, F. > soil moisture deficiency 

i 

This is the general case of a big hydrologic storm. There is 
additional streamflow due to contributions from all four components of 
runoff; i.e. channel precipitation, surface runoff, interflow, and ground- 
water flow. The hydrograph for this case is shown in Figure 1. 

2. Separation of Hydrograph Components 

The problem for a given hydrograph is how to separate the 
various components of stream flow and particularly to determine how much is 
due to surface runoff. To simplify the problem in practice, the contribu¬ 
tions due to channel precipitation and interflow are included in surface 
runoff as a single item and referred to as direct or storm runoff. This 
procedure is somewhat unfortunate for estimating contamination, but acceptable 
for the design case, since almost all of the storm runoff will be surface 
runoff anyway. Figure 1 is a typical big storm hydrograph. In it the various 
contributions to runoff are indicated: (a) surface runoff; (b) interflow; 

(c) ground-water flow; and (d) channel precipitation. The parts of a 
hydrograph are also indicated, the rising limb or concentration curve, the 
crest segment and the recession or falling limb. The time base, T, of a 



Figure 1, Hydrograph for a Large Storm, Shi 
the Various Components (after Del 


Recession 



also Flaw Contributions of 
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hydrograph is obtained by drawing a horizontal line from the point A, where 
the rising limb begins to its point of intersection B with the recession 
curve. This horizontal line may also be considered, as a first approximation, 
to be the boundary between base flow (ground water flow) and storm runoff. 

There are more sophisticated methods (l7)(2l) to separate these 
components of runoff but they consist mainly of various methods of curve 
fitting. 


Shape of the Hycrograph 


The shape of the storm hydrograph, in particular the rising 
limb, crest and early recession, is determined essentially by surface runoff. 

It is a function of sue,; storm characteristics as duration, areal distribution, 
'intensity variation of rainfall and also the shape of the basin. 


The time base of the hydrograph, T, is the time from the beginning 
of the rising limb to the time on the hydrograph when direct (or storm) 
runoff is practically zero. According to Sherman (22) the time bases of 
all hydrographs for a given basin and for storms of the same duration are 
equal. This assertion formed the basis for the unit hydrograph. The time 
base, T, may be defined by: 




+ t 

c 


( 8 ) 


where: t = duration of the storm 

K 

t = time of concentration for the drainage basin (the time 
c required for water to travel from the farthest point of 
the basin to the outlet point.) 

The effect of non-uniform areal distribution of rainfall is marked 
by a flat slope of the rising limb if most of the rainfall occurs in the 
region most remote from the basin outlet, and by a steep, rapid rise of 
the concentration curve when most of the precipitation occurs near the 
outlet. If the intensity of rainfall varies there may be more than one 

peak in the resulting hydrograph. The shape of the basin will affect the 

shape of the hydrograph; i.e. if most of the area is concentrated near the 
outlet, the peak of the hydrograph will be sooner (rapid rise), if most of .• 

the area is away from the outlet, the peak will be later (slow rising hydrograph). ; 

4. Rainfall - Runoff Relations 

The relation of runoff to rainfall is affected by so great a number 
of factors that it is obvious that no single relationship can be established 
to predict the runoff from a storm of given rainfall and duration (17). 

Nevertheless, for many years the rational approach to runoff was the use of . 

dimensionless runoff coefficients to estimate runoff as a certain percentage $ 

of rainfall. In these methods allowance was even made to subtract from the | 

data that portion of the rainfall below which there was no runoff. Other 
attempts have employed the metho ^ least squares to fit quadratic equations 
to relate rainfall to runoff, in r rational approach has some merit in the 
analysis of impervious areas such t • parking lots and airport pavements. 
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However, for natural soil conditions, it seems logical to group the data 
according to the condition of the soil before the storm ana to establish 
curves of runoff versus rainfall for each soil condition. 

5. Initial Moisture Conditions 


Several indicators of initial moisture conditions have been used: 
initial ground-water flow, soil moisture deficiency, pan evaporation data, 
total precipitation before storm. Employing the concept that soil moisture 
should decrease logarithmically with time during periods of no rain, Kohler 
and Linsley (23) devised the antecedent-precipitation index, API, which attaches 
a numerical value to the moisture conditions of the soil before a storm. 

They proposed the equation: 


I 


t 



(9) 


where: I = the initial value of API (inches of water) 

1^ = the reduced value t days later (inches of water), and 
k = a recession constant, varying from 0.85 to O .98 
6. Unit Hydrograph 

The concept of the unit hydrograph for a storm of a given duration 
for a specific basin was introduced by Sherman (22). It is based on the 
fact that the time of concentration, t , is a constant for a basin and 
therefore the time base, T = t + t. , Ss the same for all storms of the 
same duration for a given basiS. Tne shape of the unit graph is obtained 
by averaging the characteristics of many storms (on a specific basin) 
of the same duration. The area under the unit graph represents one inch 
of direct runoff from the basin. The hydrograph for any storm of the 
same duration is obtained from the unit graph by multiplying the ordinates 
of the unit graph by the storm runoff. Therefore, the unit hydrograph is 
also useful to predict peak flows. 


The above method is for storms of fairly uniform rainfall intensity 
and is limited to one-peak hydrographs. In the case of a complex storm 
with several peaks it may be possible to divide the storm into a number of 
individual hydrographs, each of which is obtained from a unit graph (17). 


Attempts have also been made at constructing the unit graph for 
storms of one duration from the graph for storms of a different duration. 
One such case is the following: Consider a storm of duration 2t^ to be 
made up of two storms each of duration tp. Shift one graph for this storm 
(this is the known unit graph) by tj* along the time base relative to the 
other and add ordinates. The resulting graph is for two inches of runoff, 
therefore, dividing these ordinates by two gives a unit graph for a storm 
of duration 2tpj. 

In the absence of streamflcw data for a particular basin, a method 
has been used by Snyder (2k) and others (17)(25) to construct synthetic 
unii graphs for the basin. 
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7. Frequency Analysis of Runoff Data 

Another basic approach that has been used for the prediction of 
runoff is based on statistical analysis of runoff data. The assumption 
is that a discharge of Q, (cfs) would occur on the average once every N 
years. If it were decided that, say, a 100-year event would be sufficiently 
probable to take into consideration, the runoff corresponding to this 
probability would be the amount chosen to work with. 

8. Rainfall Intensity - Duration - Frequency Curves 

Since there are not many locations that have runoff data for 
periods long enough to avoid extensive extrapolation in the prediction 
of even a 100-yr. runoff event, appeal has sometimes been made to the 
usually much longer rainfall intensity - duration data available. This 
method assumes essentially that a 100-yr. rainfall event produces a 100-yr. 
runoff event. The previous discussion of the relation between rainfall and 
runoff should make clear the dangers inherent in this assumption. 

F. Drainage Basin 

The drainage basin or watershed, of a stream is that area that 
contributes runoff to the stream; it is the natural unit for many 
hydrological studies. The boundaries of a basin are divides or ridges 
which separate it from adjacent basins. A basin has a single outflow 
point, either where the stream cuts through a divide or where the 
stream reaches the ocean. 

First of all, it is necessary to consider the physical characteristics 
of the basin itself or its morphology (26)(27). Some of the relevant 
characteristics are: (a) drainage area of the basin (the area contributing 
to runoff), A^, and (b) drainage density (stream length in the basin, 
area) D^ = li/A^, where £l = total stream length in the basin. Values 
of D,j vary from less than 1 mile/sq. mile for a poorly-drained basin to 
about 5 mil/sq, mile for a well-drained basin. 

Streams may be typed by their flow constancy as follows: (a) perennial 
(at all ti.ies above surface flow except during time of extreme drought), 

(b) intermittent (mostly above surface but sometimes dry due to evapora¬ 
tion and bank storage), and (c) ephemeral (from flash runoff or snowmelt 
only). 

To compute the time distribution of runoff the distribution of area 
with distance from the outflow station is considered. On a basin map 
isopleths (lines of equal distance along streams from the outflow point of 
the basin) are drawn and the percentages of the drainage area beyond an 
isopleth versus the isopleth value in miles are plotted. This gives the 
distribution of drainage area as a function of distance from the outflow 
point. 

Another distribution curve, a plot of elevation vs. area, may be 
constructed from a topographic map and is useful in comparing drainage 
basins. The curve is a plot of the average elevations vs. the areas 
in percentage of the total area. Either the median elevation, such 
that one-half of the area is above and one-half below, or the mean 
elevation, which is a weighted average, may be used. 
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The stream profile shows the elevation of the main bed of a stream as 
a function of its distance from the basin outflow point. The gross slope 
of a stream between any two points is the total fall between the points 
divided by the stream length. The mean slope is constructed by drawing 
a straight line such that the areas enclosed above and below the stream 
profile are equal. The stream profile and the mean slope are specific 
characteristics of the stream channel and are not suitable for use as 
parameters describing the slope of the overall basin. The average land 
slope of a basin may be determined by a method of superimposing grids 
^ver a topographic map according to Horton (28). 
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V. ION EXCHANGE IN SOILS 


The interaction of fallout particles deposited on the watershed with 
the soil is of great importance in this study, lfce deposited nuclides 
that are dissolved by surface runoff may undergo ion exchange with the 
surface layers of the soil. Certain of the isotopes are held quite tightly 
by soils, such as cesium by clays and hence are not expected to appear in 
runoff. Others, such as strontium are generally not sorbed very strongly. 
Reports of relevant experimental work from the literature, are summarized 
in Section A of this Chapter. 

A simplified mathematical model of uptake by soils is presented in 
Section B. This model has the virtue of providing an analytic solution. 

It is intended as a first approximation since it was derived for the 
case of isotopic exchange. Strictly speaking this case obtains rarely 
in practice. Equations incorporating more general ion exchange processes 
should be used. These equations entail much greater mathematical complexities 
and the analysis is deferred to a subsequent report. 

It should further be noted that this model is restricted to a one¬ 
dimensional case and to moisture-saturated soil. In practice, the soil 
is not saturated during at least part of the period of infiltration. 

This model or a subsequent one, is to be incorporated into an analysis 
that will provide the overall picture of transport of radionuclides in 
overland flow and infiltration and will lead to an accurate description 
of the distribution of the activity between these two categories. 

A. Literature Survey of Experimental Results 

Many nuclides have a strong affinity for cation exchangers, such as 
clays and soils, and may be retained by the soil rather than removed by 
surface runoff. The amount actually removed by runoff depends on the 
chemical composition of the soil and the amount of organic matter it 
contains (29)(30)(31)(32)(33)(3*0 • Reports in the literature generally 
agree that pH greatly affects nuclide sorption, with the best sorption 
occuring at the higher pH values, as illustrated by Table II for 
Savannah River Plant Soil, according to Prout (35): 


TABLE II 

Exchange Capacity of Soil (Savannah River Sediment) 


£2 


meq/100 g Soil 


meq/100 g Clay 


4 

5 

6 

7 

8 

9 

10 


0,012 

0.108 

0.372 

0.504 

0.600 

0.744 

1.24 


0.06 

0.54 

1.86 

2.52 

3.00 

3.72 

6.18 
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Results of studies conducted with the more common clay minerals 
show that they have different affinities for different nuclides (29)(35) 
( 36 -te). As presented by Lacy (43), most clays have preference for competing 
ions in the following order: 


Ba ++ > Sr** > Ca** > Mg ++ 
and 

Cs + > Rb + > NHj, + > K + > Na + > Li + 

Additional data on the cation exchange capacity of a number of clays (44) 
has been tabulated in Table III below. 


TABLE III 

Important Properties of Some Clay Minerals 


Group and Lattice Type 

Name 

Cation exchange 
capacity 
(meq/100 g ) 

Approximate 
size range 
(**). _ 

KAOLINITE group 
(1:1 lattice) 

kaolinite 

3-15 

0.5-2 


balloysite 

1*0-50 

0.04-0.2 

M0N3M0RILL0NI IE group 

(2:1 expanding lattice) 

raontmorillonite 

80-125 

0.01-0.1 


beidellite 

60-90 

0.05-0.5 


nontronite 

60-70 

0.01-0.1 


saponite 

20-30 

0.01-0.1 

ILLITE or HYDRATED 

MICA group (2:1 non¬ 
expanding lattice) 

illita 

20-1*0 



pyrophyllite 

- 

- 


As the result of the initial literature survey, experimental data for 
specific cases have been summarized in Table IV, as shown on pages 25 and 26. 
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TABLE IV 

Stannary of Experimental Work on Ion Exchange Uptake in Soils (Sr and Cs 


Element 

Soil 

Results 

Reference 

Sr-90 

(Los Alamos) tuff 

Hot tightly bcarnd; easily 
replaceable 

45 

Sr 

Savannah River Plant 
Soil 

High H strongly inhibits 
adsorption} Ha 4 and 
inhibit adsorption above a 
pH of 8; Ca 24 seriously 
inhibits adsorption; Ha and 
A1 compete with Sr for ex¬ 
change sites but not as much 
as with Cs 

35 

Sr 

Vermiculite 

(P0h) J improved sorption of 
Sr when added in concentra¬ 
tions up to 100 ppm 

36 

Sr 

Montmorillonite 

Sorption reduced by the 

36 


following ions in this order. 

Ba 2+ > Cs2+> Mg2 + > h+> nh£> 

K 4 ^- Na 4 . Indicating that Mont¬ 
morillonite sorbs Sr by ion 
exchange 

2 + 

Sr Montmorillonite, Ca increases sorption up 36 

vermiculite, varis- to a concentration of lOOOppm 

cite, Tenn. rock 
phosphate, Fla. pebble 
phosphate 

Sr Coalinga, Dominguez, H + concentration had a marked *+0 

Lost Hills-Asbestos, effect e.g. reducing pH caused 

Richfield, Rosecrans, significant reductions in Sr 

Halloysite Clay-Asbes- uptake. Equilibrium constants 

tos, Yolo Soils-Asbestos were determined for Sr-Ca which 

was around 1.3. It was depend¬ 
ent on the H 4 concentration. 

Sr Rosecrans sand With Na in the liquid phase 

Sr-Ca equilibrium constant 
was 1.79 

Sr Kaolinite, halloysite, The concentration of Ca de- 46 

montmorillonite, creased Sr sorption more than 

illite did Na concentration. Author 

concludes that to predict sorp¬ 
tion of Sr by Clays it is 
necessary to consider Na-Ca 
ratio in the water as well as 
total cation concentration 
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Summary 

Element 

Sr-90 


Sr-90 

Cs-137 

Cs 

Cs 


Cs-137 


# 


TABLE IV (Cont'd) 


Experimental Work on Ion Exchange Uptake in Soils (Sr and Cs) 


Soil Results Reference 

World-Wide soil dis- There was little downward 47 

tribution,geographically movement of the Sr-90 through 
indexed soils. It takes several years 

of heavy rainfall for Sr-90 
to be moved even a few inches 
downward 

Clay 97$ of the Sr was removed by 48 

50 ml of leachate when leached 
by 100 ml of HgO and NK^ acetate 

(Los Alamos) tuff Tightly bound 45 


Savannah River Plant 
Soil 


Clays 


Clay 


Decrease in adsorption at 35 

pH below 6 due to H*; high 
pH decrease in adsorption 
due to Na' 1 ' 

As Ca or Na concentration is 46 
increased,Cs sorption is 
decreased. The author concludes 
that as a first approximation 
the sorption of small quantities 
of Cs by clays would be depend¬ 
ent only on total cation con¬ 
centration in the water 


Only 51$ removed in the first 48 
50ml of leachate when leached 
by 100 ml of HgO and HFfy acetate 

total activity of Cs-137 leached » 

was that of Sr-90 j 

I 


i 
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Model for Uptake of Radioactivity by Runoff - 
Diffusion Into a Rod from a Solution of Finite Volume 


The model, except for the finite volume of the solution, is quite 
analogous to sorption by clays. In the experiment, described by equation (ID) 
below, a rod, initially free of radioactivity, is insulated except for one 
end which is exposed to a radioactive solution circulating rapidly past 
the rod. The decrease in activity of the solution Q/Qq, is given by: 


SJil, / 


erfc (f) 


where: 



2 

A = cross-sectional area of rod (cm ) 

O 

V = volume of solution (cnr) 


C* -• equilibrium concentration of radioactive species 
in solution (A22JL) 
cm3 

C* = equilibrium concentration of radioactive species 
in the rod ^ ons ) 
cm3 

It should be noted that in the model the amount of activity in the 
solution is initially fixed, while in the case of uptake by soil, the 
solid is continuously exposed to a constant concentration of radioactive 
solution. 


Clearly the above model differs from the actual situation in uptake 
by soil from runoff in many aspects. The model does have the virtue of 
having an analytic solution. Tne authors intend this model only as a 
first approximation} the actual equations and boundary conditions for 
ion exchange (as distinct from isotopic exchange as considered in the 
model) should be solved numerically on a computer. 


Hie strong analogy between the present problem and the operation of an 
ion-exchange column has also been noted. 
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VI. PLANT UPTAKE 


A. General 

In their course of travel before reaching a water source, fallout nuclides 
may also be taken up by the plant cover of the watershed. It is necessary 
to know the behavior of the various nuclides with respect to plant uptake, 
as this will affect the final amount of activity reaching water supplies. 

Fallout nuclides may be taken up by plants in two ways: 

(a) The may land directly on plant parts and be taken in through 
the leaves, or,(b) they may be taken up by the roots. 

Nuclides which land on the plants may also be washed off the plant 
by rain and enter the water system by runoff. The uptake by plants, 
however, is important to estimate the total amount available for future 
runoff (37)(^+l)(^2)(^9)* The composition of the soil has a distinct 
effect on plant uptake of nuclides (37)(38)(39)(50)(51)(52)(53)(5*0(55). 

B. Uptake of Sr-89, 90 

Schulz, et al.(51) ran culture studies on plant uptake of Sr-90. 

They found greater uptake by the plant than extraction by water. They 
stated that conditions in the soil may be different than the culture 
and the amounts of Sr extracted may not be the amounts available to plants. 
Schulz, et al., also indicate that some reports have found Sr to be in 
unexchangeable forms which could be due to the type of soil and length 
of time Sr has been in the ground. 

Libby (37) found that plants contain about twice the specific Sr-90 
relative to Ca in soils on which they grew. This may be due to fallout 
on plant surfaces. Menzel, et al., (56) found that the Sr-90/Ca ratio 
of cowpeas was "approximately inversely proportional to the available Ca 
in the soil ..." 

Various investigators have found that addition of Ca to soils as 
gypsum, or lime, decreased the Sr-89, 90 uptake, except in soils already 
calcareous (55). Nishita, et al. (30)(3l), using various organic fertilizing 
materials, ran tests on Sr-90 uptake. Lettuce had the greatest effect with 
10 gm lettuce/100 gm soil reducing the Sr-90 content by 20 to 40$ of plants 
grown on soil with no addition. This decrease due to the addition of organic 
materials may be explained in several ways such as: Microbial activity, 
changes in composition and amount of soil air, decreasing the ratio of 
divalent to nonvalent cations, and the amounts of mineral nutrients being 
increased (55). 

C. Cs-137 Uptake 

Nishita, et al. (39), found the Cs-137 uptake was dependent on the 
potassium concentration. The greater the potassium concentration in the 
soil the less Cs-137 was taken up. The addition of small amounts of non- 
radioactive Cs-137 greatly increased the uptake of Cs-137. The ratio of 


j 


i 
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Cs/K in uptake depended on the soil type. Handley, et al. (42), used a 
greater number of ions with the Cs-137 and found X, Rb, NHj^ and Cs bad 
large effects on Cs-137 with K and Rb having the greatest deterring effect. 
I'iddleton, et al. (4l), observed that barley absorbed K-42 to a greater 
extent than Cs-137 relative to the concentration of the external solution 
in the culture. Jackson, et al. (38), found that different salts had 
different effects. Some uptake results were quite different from those 
of Handley, et al. (42), who found Rb and NH^ decreased Cs-137 uptake, while 
Jackson found the opposite as shown in Table V. 

D. Comparison of Uptake 

It has been found that plant uptake is in the following order: Sr-89, 

90 » 1-131 > Ba-l40 > Cs-137 and Ru-106. The uptake of Sr-89, 90 was 
from 0.5 to 5$ of the amount in soil, while the other nuclides are taken 
up in smaller percentages (57)(58)(?9)• 

Klechkozsky (60), using Sr-90, Ce-l44 and Ru-106, observed that strontium 
uptake was reduced 20$ when the soil war limed. This reduction was pro¬ 
bably due to the fact that the solubilities of Ce and Ru were decreased 
with an increase of pH. The presence of NH^NO^ increased the uptake of 
all cations. 

In several Russian experiments using Sr-89, 90, Ce-l44 and Ru-106 it 
was found that uptake was decreased in all instances when lime and organic 
matter were added. Also, legumes take up the nuclides to a greater extent 
than do cereals. Strontium is taken up in the greatest amount by the plants 
and it is most radically reduced by the additives. In one experiment, 
when only organic matter was added, the decrease was greater than when 
only lime was added, but less than when both lime and organic matter were 

used (32)(33). 

To determine if an increase of the sorptive capacity of the soil is 
caused by addition of lime and organic matter to the soil, Guliakin, et al., 
(32)(33) performed experiments in which one gram of soil was placed in a 
test tube to which 10 ml of a solution of measured radioactivity was added. 

The soil was shaken and centrifuged and the activity in the solution 
measured again. This was done again with salts, lime and organic matter 
added at various points of the experiment and in various combinations. 

The results demonstrate that lime and organic matter had no marked 
effect on the amount of Sr, Ce, and Ru sorbed by the soil. There is a 
reduction in the displacement of Sr-89, 90 and Ce-l44 from the absorbed 
state when lime or organic matter is introduced into the soil even by calcium 
and potassium salt. The salts displaced almost no absorbed Ru-106 even 
when lime and organic matter were introduced into the soil. Ru-106 was 
reduced in plants whose soil had been treated with lime and organic matter. 

It may be concluded from these observations that a plant has greater 
sensitivity to a change in nuclide sorption than shown by these experiments (32). 

Another experiment seemed tc indicate that plant uptake is greater in 
sandy loam than in loamy soil. This phenomenon is usually due to the 
availability and greater abundance of plant nutrients in loamy’ soil. The 
reduction of uptake when lime and organic matter is added is greater for 
the sandy loam than for loamy soil (33)* 

As a part of these studies, a critical review of recent literature 
on the chemical and exchange properties of fallout nuclides in various 
media was carried out and over 150 references abstracted. 
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TABLE V 

Ion Effects on Nuclide Uptake by Plants 
According to Jackson (38) 


(*) 

Salts Added' ' 

Water Extractable Cs-137 
cpm (x 10~3)/pot of soil 

Uptake of Cs-137 
cpm (x 10“3)/gm of plant tissue 

none 

0.70 

± 0.10 

24.76 

- 2.50 

KH^Cl 

13.72 

± 2.05 

44.69 

± 0.98 

nh 4 no 3 

9.93 

i 1.71 

4l.4o 

i 4.53 

(m 4 ) a so It 

11.24 

± 1.71 

44.15 

± 1.00 

KC1 

5.09 

± 1.01 

1.23 

-.0.15 

“°3 

4.19 

± 0.52 

2.00 

± 0.62 

w 

4.91 

i 0.63 

1.19 

± 0.08 

NaCl 

1.08 

- 0.30 

33.84 

± 3.09 

CaC^ 

1.64 

- 0.19 

34.68 

- 3.73 

MgClg 

1.80 

± 0.20 

34.94 

i 1.58 


^ ^The salts were added at 1 meq/100 gm of soil 


I 

& 

I 
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VI. STATUS OF COMPUTER PROGRAMS 


The development of computer programs has followed along the lines 
suggested in the scope of this research project. An estimation of tne 
number of atoms per unit area, N(A), is a modification of the evaluation 
of fallout intensity, 1(1), at a given location. To this end, a 
program for interpolation of nuclide solubility contour ratios from data 
supplied by Miller (6l) has been developed. This program is independently 
useful in that it provides the value of the contour ratio for each downwind 
distance, or particle size parameter, a. Furthermore, the incorporation 
of this program into a MAIN program such that the final output will be 
in atoms per area over a dcwnwind-crosswind surface is nearly completed. 

A computer flow diagram appears later in this section of the report. 

The evaluation of atoms per area for the cumulative effects of a 
poly-weapon situtation represents an extension of the single weapon atom 
concentration program. Such a program is necessary for the evaluation 
of contamination from large-scale attacks affecting watersheds and reservoirs. 
The program previously developed for multiple weapon evaluation of intensity 
will represent only a minor portion of the expanded version. A new multiple 
yield program to compute the total atom concentration at any point or for 
the entire watershed is now in the developmental stages. 

A. Development of Computer Program to Determine Atom Concentration, N(A), 

Over Any Surface 

In the program entitled Soluble Nuclide Contour Ratios, previously 
presented (l), the computation of N/l (or N^i(l)) was expected ultimately. 

In the meantime these values have been computed by Miller. Therefore, to 
facilitate the evaluation of contaminated water supplies, it has been decided 
to use the nuclide solubility contour ratio data provided by Miller (6l). 

For the evaluation of water contamination due to a surface blast of 5, 10, 
or 20 MT a tedious procedure had been followed. The difficulty was due to 
the lack of a completely successful computer program to evaluate the 
nuclide solubility contour ratios. By incorporating the contour ratios 
(N/I) computed by Miller into the existing program (62) for intensity 
evaluation (i), the number of soluble atoms (N) can be found readily. 

The previously developed program for estimating fallout intensity at any 
location has been streamlined and modified by incorporating the soluble 
nuclide contour ratio, N$i(l), into the MAIN program. Modifications 
include a subroutine to interpolate for the desired values of N2i(l), 
equivalent to N/l, and a multiplication of these contour ratios by the 
computed intensity values, or 


N(A)-= I(l)!^ i (l) = I x | = N 

With this modification the final output will be in terms of atoms/area. 
The outline of the entire program, including this modification, may be 
seen in the flow diagram, Figure 2. 
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The table from Miller (6l) supplies only values of N?i(l) for 1, li), 

100 MT weapons. For this reason it was necessary to interpolate for 
5 and 20 MT yields. The Vf (converted to a*) values in the table are not 
sufficient for an accurate determination of contamination over large water¬ 
shed areas. In the flew diagram given below, the SUBROUTINE INTER? is a 
method of linear interpolation for the necessary a values. The subroutine 
ha 3 been tested separately in a slightly different form and is presented 
with its flow diagram later in this report as Figures 3 and 4. 

The use of a computer program for the direct evaluation of N '(atoms/area) 
has two justifications; 

(1) The large number of coordinate points within a given watershed 
area may be quite large. Therefore,fche number of calculations for intensity 
and the interpolations for nuclide contour ratio values would be quite 
voluminous. 

(2) By proper selection of coordinates, the computation may be made 
as accurate as initial assumptions allow. The desired degree of accuracy 
is achieved by decreasing the square area of evaluation, 

In the following pages, the flow diagram for the entire computer program, 
that for the interpolation subroutine and finally the actual computer 
program with some output are presented. 
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READ INPUT 

1 ) 

2 ) 

3) 

4 ) 

5) 

6 ) 



Over Any Surface 

Note: a* are values corresponding to Vf selected by Miller (6l). ,, 

a are values calculated according to the formula, c#=*x/hsXi/(l.68xl0 4 Vr . 


weapon yield (W) 
downwind distances (X) 
contour profile points 
crosswind distances (Y) 
oc* according to Miller 
1^(1) according to Miller 
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Figure 3. Flow Diagram for the Interpolation of Nuclide Solubility 

Contour Ratio, N$i(l), for Values of a between 110 and 0.667' / 


4 


{*) Note: corresponding to values of V f from 0.2 to 33.0 














C INTERPOLATION OF NUCLIDE SOLUBILITY CONTOUR RATIOS 

DIMENSION A < 8 ) * ANFW(IOO). SR89!8»3)» SR90(6»3)* RU106C8*3)* 

1 XI131(8*3)* CS137 ! 8 * 3 ) * BA140!8»3)* U < 3 )* V C 31 * W(3>* X C 3)* 

2 Y ( 3 ) ♦ 2(3) 

RFAD 1* N»M»{A(I)*I=1*N)#(ANEW(tO *K=1*M) 

READ 2 * ((SR89{I * J)»J=1.3)*<SR90{I»J)♦J=1»3)»IRU106(I *J)*J*1 *3 ) ♦ 

1(XI131CI*J)»J=1»3)» <CS137(I♦J)♦J*1» 3)•(BA140(I♦ J)*J*1*3)* I*1»N) 

DO 11 K= 1 »M 
DO 10 1 = 1,N 

IF(ANEW(K)-A CI)) 10,13*12 

10 CONTINUE 

13 PRINT 21* ANEW(K) 

21 FORMAT(22H0ANEW IS A DATA POINT , F10.5) 

GO TO 11 

12 CALL INTER!I.A.ANEW,SR89»SR90.RU106*X 1131*CS137,BA140* 
1U,V.W,X*Y,Z) 

PRINT 20, ANEW(K) ♦ U* V* W« X, Y# 2 

11 CONTINUE 

20 FORMAT!40H1THF ALPHA W' ARE INTERPOLATING ABOUT IS* F10.5// 

125H0THF VALUES FOR SR89 ARE* 3F10.3//25HOTHE VALUES FOR SR90 AR 
2E * 3F10.3//26H0THE VALUES FOR RU106 ARE* 3F10.3// 26H0THE VAL 

3UES FOR XI 131 ARE, 3F10.3//26H0THE VALUES FOR CS137 ARE* 3F10 
4.3//26H0THF VALUES FOR BA140 ARE. 3F10.3//////) 

1 r .ORMAT( 2I5/(8F10.4) ) 

2 FORMAT(6F12.4) 

CALL EXIT 

END 


SUBROUTINE INTER!I,K,A♦ANEW,SR89»SR90,RU106*X1131 *CS137,BA140,U. 
U ♦ W »X ♦ Y , l ) 

DIMENSION a c ft j * A!L W ( 1 0 0), SR89<8*3)» SR90!8»3J* RU106(8»3)» 
Xi:p .3), C c 3 7 ( 8 ♦^), PA 140{R, 3 ), Ui3), V(3)» W(3)* X!3)♦ 

? .' ( 3 ) ♦ Z ( 3 ) 

3=(ANEW(<)-A( I) ) 

C*(a(1-15-A! I ) ) 

DO 10 J=1,3 

10 U(J) = (SR89(I-1»J)-SR89(I*J))*B/C+SR89< I*J) 

DO 11 J= 1,3 

11 V!J) = (SR90{I-1*J!-SP90(I»J))#B/C+SR90(I *J) 

DO 12 J=1*3 

1? W(J)=(RU106(1-1.J5-RU106!I♦J))*B/C+RU106!I»J) 

DO 13 J=1,3 

13 X(J) = (XI131(I-1*J)“XP31 (I ♦ J) )*B/C+XI 131(1 ,J) 

DO 14 J=1» 3 

14 Y(j:=(CS137(I-1*J)-CS "<I«J))*B/C+CS137(I*J) 

DO 15 J=1» 3 

15 It J)=(BA140(I-1,J)-BA140(I*J))*B/C+BA140<I♦J) 

RETURN 

END 

* DATA 


Figure k, Computer Program for Interpolation of Nuclide Solubility Contour 
Ratio, n£i( 1), for Values of a between 110 and O .667 (#) 

(*) Note: corresponding to values of V f from 0.2 to 33.0 
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5>'T 

1GM i 

20M i 

THE 

values 

FOR 

SRfi9 ARE 

0*470 

0.023 

0 • 0 2 u 

THE 

VALUES 

FOR 

SR90 ARE 

0*380 

0.197 

0.145 

THE 

values 

FOR 

RU10A ARE 

0.125 

0.095 

0.0^5 

THE 

values 

FOR 

XI131 ARE 

0.08 5 

0.055 

U.006 

TEE 

values 

FOR 

CS137 ARE 

0.100 

0.075 

0.040 

THE 

values 

FOR 

bA14C ARt 

1.165 

0.584 

0.503 


THE ALPHA w£ ARE INTERPOLATING AbO-Ji IS 68.75 






5MT 

10MT 

20 N'T 

THE 

values 

FOR 

SRA9 ARE 

4.1 AC 

4 . 5 A 0 

4.5K1 

THE 

values 

EOR 

SR90 ARE 

0.000 

*.'>4 5 

5.099 

THE 

values 

FOR 

RU106 ARE 

2.220 

2*360 

2.357 

THE 

values 

FOR 

X11 31 ARE 

u.923 

1.015 

1.u54 

THE 

values 

FOR 

CS13: ARE 

o« obo 

7.630 

/• 

THE 

values 

FOR 

om140 mRc 

6.495 

5.705 

0 . 3 ! o 


Figure 4. (cont'd) Computer Program for Interpolation of Nuclide 

Solubility Contour Ratio, n2j( 1), for Values 
of a between 110 and 0.667 (*) 

(*) Note: corresponding to values of V from 
0.2 to 33.0 f 
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3. Development of Multiple Weapon Program 

Incorporation of contour ratio data into the intensity evaluation 
program implies a similar extension of the multiple weapon intensity 
program. The multiple weapon program presented in the previous report (1), 
was inefficient and required considerable streamlining. The new multiple 
weapon program necessarily follows a successful single weapon yield 
program. Though not completed, the program is expected to follow the 
guidelines of the following flow diagram. 

It sr.ould be noted that Nri(l) values are interpolated in a sub¬ 
routine from the data supplied by Miller (6l). Also, there are two 
subprograms available for the translation of axes. One of these uses 
the rotation of axes presented in the previous report (1). The ether 
subroutine assumes parallel axes. In some instances a method of 
parallel translation may be more appropriate. Since ,?et streams ana 
upper air wind patterns do not change radically, two target cities may 
well experience the same wind direction. In fact, if two differing ground 
zero blasts are cc have a noticeable effect on a watershed area, they 
will certainly be within a six Hundred mile range of each ether. It 
can be expected that in such e case tne wind directions at both sites would 
be approximately equal. In this case the parallel translation of axes 
would be the more appropriate method of analysis. 


1 
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START 


READ INPUT 


1) downwind distances (X) 

2) crosswind distances (Y) 

3) cc* according to Miller 

4) N^i(l) according to Miller 


SELECT the data 
for next weapon 
considered 


READ INPUT 

1) weapon yield (W) 

2 ) contour profile points 

3 ) angle of rotation 

4) perpendicular distances 
between axes 


TRANSFER NECESSARY 
VALUES: 

1) weapon yield - 

2 ) downwind distances 

3 ) profile points 

4) crosswind distances 


COMPUTATION of arithmetic 
expressions which reoccur in 
the evaluation of l(l) 


CALL for SUBROITINE INTENS 


CALL for SUBROUTINE INTERP 


EVALUATION of Ifr^) 
for i = 1, N and j ~ < ’1,M 


TRANSFER NECESSARY 
VALUES: 

1) a* according to Miller 

2 ) N$i(l) according to 

Miller 

3 ) downwind distance 


RETURN 


TRANSFER NECESSARY 

VALUES: 

1) downwind distances 

2 ) crosswind distances 

3 ) other parameters 
depend on which 
method of translation 
is used 


TRANSLATE each coordi¬ 
nate point by the 
appropriate formula for 
the method used 


CALL for SUBROUTINE TRANS 
(Translation) 


MULTIPLY Nj^l) by l(l) 
for (Xj_,Yj, isotope) and 
sum with previous N(A) 


PRINT 

atom concentration, N(A) 


COMPUTATION ofa's 
a = X i /(l.68xl0 5 +W 0 ' 164 ) 



DETERMINE between which 
: successive o*'s the 
| computed a lies 


INTERPOLATE the N?i(l) 
value for the six 
isotopes 


RETURN 


RETURN 


Figure 5. Flow Diagram for Atom Concentration at Any Location for Multiple 
Weapons Computer Program 
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VIII. PRELIMINARY EVALUATION OF CONTAMINATION OF THE PROVIDENCE, 
R.I., WATER SUPPLY 


A. Selection of Attack and Orientation of Fallout Model 


A preliminary evaluation of fallout contamination in the Providence, 

Rhode Island, water supply system was carried out during the initial 
phase of this research program in preparation for the Five-City Study, 
sponsored by the Office of Civil Defense. 

In consideration of the limitations in time and resources for this 
phase of the study and in the absence of an operational multiple weapon 
program, it was decided to follow the previously reported Technical 
Operations, Inc. attack model (63) to provide the fallout contamination over 
the Providence, R.I. watershed, with an assigned wind velocity of 15 mph. 
Ultimately the study will examine e variety of attack conditions. It 
followed logically zo select Springfield, Massachusetts, as the target 
city because of the predominant wind direction and to supply an additional 
10 MT weapon to approximate the effect of the other weapons in the attack 
model. Location of the weapons with ground zero at Springfield, 

Massachusetts, probably accounts for over 80$ of the fallout on the Providence, 
R.I., watershed. Therefore, the fallout pattern downwind axis was 
oriented in an east-southeast direction, as illustrated in Figure 6, 
for the following reasons: 

1) The predominant wind direction at Providence, R.I., is west- 
northwest. Since the fallout level at any particular point downwind 
is very sensitive to wind direction and the entire watershed is 
located due west of the city, location of ground zero at Providence 
would not provide fallout contamination of the city's system of 
water supply reservoirs. 

2) The relevant geographic industrial and military target areas 
presented in the attack model (63) and consideration of their 
locations relative to the Providence watershed led to the 
selection of cne 10 MT ar'’ 20 MT weapon detonated at Springfield, 
Massachusetts, situated 60 miles west-northwest or Providence, 

Rhode Island. 

It should be pointed out that these estimates of water contamination, 
for six biologically important radionuclides, as well as those reported 
for other target cities previously (l), consider the effect from an attack 
on a single target area only. In any general, multiple-weapon attack on the 
continental United States, fallout from weapons detonated over other target 
areas would result in overlapping of coverage and increase the contamination 
level over reservoirs and in water supplies. An unclassified attack pattern 
for use in the initial analysis of the Five-City Study has been evolved, 
but this Guide (61) was not received until after the work covered in this 
report had been completed. 




- 4o - 



Figure 6. Location of Ground Zero for Evaluation of Providence Water 
•Supply Contamination 




B. Providence Water Supply System 


The present system of water supply for the City of Providence was 
constructed under the supervision of a commission of seven members, known 
as the Water Supply Board, created in accordance with Chapter 1278 of 
the Public Laws of Rhode Island and approved on April 21, 1915* The 
Cities of Providence and Cranston along with the Towns of Johnston and 
North Providence are supplied through a distribution system owned and 
maintained by Providence. Parts of the City of Warwick and -portions of the 
Towns of West Warwick and Coventry served by the Kent County Water Authority, 
along with the East Smithfield Water District, are furnished water on a 
wholesale basis but own the distribution system within their respective 
communities. Other towns entitled to Providence Water under legislative 
acts, but not being supplied at present, are Scituate, Foster, and 
Glocester. These cities and towns represent about 377 sq. mi., or about 3 6 $ 
of the land area of the State of Rhode Island. 

1. Water Consumption 

Based on the i 960 census, the Providence Water Supply Board 
provides water to 45$ of the population of the state of Rhode Island. 

In 1962 , the eight cities and towns receiving Providence Water used 
45.72 million gallons per day. Coincidental with the severe drought, the 
average daily water consumption increased to an all-time record for the 
water year ending September 30, 19^5. During this period, the cities and 
towns supplied from the Providence system consumed 54.6 million gallons 
per day(*). According to the 1965 Annual Report of the Providence Water 
Supply Board ( 65 ) a total population of approximately 383,635 is served by the 
system. The rate of daily water consumption has increased from approximately 
119 gallons per capita in 1962 \o 139 gallons per capita in 1965 , or 
almost 17 $, which may be at least partially attributable to the prolonged 
drought experienced in the Northeast. 

2. Watershed 


Providence obtains its water from a surface supply located on the 
north branch of the Pawtuxet River. The total watershed area covers 9^*8 sq.mi. 
as shown in Figure 7 . The watershed area is about five times the area of 
the City of Providence. This drainage basin represents approximately 9$ of 
the land area of the State of Rhode Island. The City owns 23.93 sq. mi., 
or slightly over 25 $ of the land in the drainage area, or about 5 sq. mi. 
more than the 18.91 sq. mi. area of the City of Providence ( 66 ). 


(*)The total annual draft from the Scituate Watershed for 1965 was 

25.06 billion gallons, or an average of 68.65 million gallons per day. 
The average daily draft for water supply purposes was 54.6 million 

t allons a’'d the difference cf 14.05 million gallons per day was 
ischargeu iitc the north branch of the Pawtuxet River. 






Figure 7. WATERSHED OF PROVIDENCE, R.I. WATER SUPPLY 
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3. Reservoir System , 

; i 

All surface water from the watershed is ultimately collected in 
the main reservoir, known as the Scituate Reservoir. Five smaller 
reservoirs: Moswansicut, Regulating, Ponaganset, Barden and Westeonneug 
are tributary to the main reservoir. Further information on these reservoirs 
is presented in Table VI below. 


TABLE VI 

Providence, Rhode Island, Water Supply Reservoirs 


Reservoir 
(Name) 

Watershed Area 
(sq. mi.) 

Storage Capacity 
(million gallons) 

Water Surface 
(sq. mi.) 

Spillway Elevation 
(feet) 

Scituate 

92.8 

37,011 

5.30 

284.01 

Moswansicut 

3.9 

1,701 

0.44 

301.90 

Regulating 

22.3 

128 

0.38 

285.50 

Ponaganset 

2.1 

724 

0.36 

633.05 

Barden 

33.0 

053 

0.38 

345.10 

Westconnaug 

4.0 

453 

0.27 

454.17 


All five small reservoirs were originally owned and controlled by 
mills located along the Pawtuxet River and wer* acquired by the City under 
the provisions of the 1915 Water Act. Gross storage in all reservoirs 
totals 41.268 billion gallons but the dead storage amounts to 1.522 billion 
gallons, which leaves a total available storage of 39.746 billion gallons. 

At the end of the water year, September 30, 19^5, the combined 
storage was 29.407 billion gallons, or 71.3$ of capacity. The maximum 
combined storage occurred on May 2, 1965, when 30.55 billion gallons, or 
93.4$ of capacity, were impounded. The total annual draft from the 
Scituate Watershed was 25.056 billion gallons, or an average of 68.65 
million gallons per day. The total annual draft for water supply purposes 
was 19.93 billion gallons, or an average of 54.6 million gallons per day. 

4. Hydrology 


The average annual rainfall on the watershed is 48.42 inches, based 
on the 47-year average (1915-1962), with a 66.28 in. yearly maximum (1950) 
and a 33.43 in. yearly minimum (1957). The average yearly runoff, or 
water actually collected in the reservoirs, based on the same period of 
record, is 25.12 inches. Every inch of runoff over the 92.0 sq. mi. 
watershed is equivalent to a volume of 1,612,750,000 gallons. Multiply¬ 
ing this figure by the long-term average runoff of 25.12 in, and dividing 
by 365 days shows an average yield of 110,990,000 gallong daily, or about 
twice the average quantity of raw water delivered daily to the Water 
Purification Works. However, the estimated safe yield of the Scituate 
Supply is 84,020,000 gallons per day. 
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Rainfall on the 92.8 sq. mi. Scituate Watershed, due to the extended 
severe drought, was only 38.13 in. for the year ended September 30 > 19^5. 

This was the fourth lowest annual rainfall experienced during the 50-year 
(I916-I965) period of record; it was 9*91 in. less than the long-term average 
of 46.04 inches. 

The runoff for the year 1965 totalled 14.26 in. which was 10.40 in. 
less than the 50-year (1916-1965) average of 24.66 inches. It was the third 
lowest annual runoff during the 50 years of record. 

For the year 1965, the yield from the Scituate Watershed was 
22.99 billion gallons, or an average of 62.99 million gallons per day, 
which was 5.65 million gallons per day less than the average daily draft, 
and 45.897 million gallons per day less than the 108.89 million gallons 
per day average yield for the 50-year period 1916 through 1965. 

5. Watershed .Management 

The City of Providence has a highly developed system of watershed 
maintenance and management. All forestry operations on the watershed are 
under the management of a professional forester who supervises the work on 
over 11,000 acres of City-owned forest land that surrounds the main Scituate 
Reservoir and the five smaller reservoirs. Practically all the arable areas 
have been planted with conifers such as White Pine, Red Pine, Scotch Pine, 
Austrian Pine, Jack Pine, White Spruce, Nonray Spruce, etc. Much of the 
wooded area, particularly that upon which low-quality oak and hardwood grew, 
has been underplanted with coniferous species. Approximately 7,000,000 
trees have been planted on watershed lands owned by the City. 

Forest cover is essential on the Scituate Reservoir watersned 
for the storage and supply of high quality water. A thick leaf litter on 
the forest floor serves as a soft, spongy surface to absorb the rain and 
melting snow. The undisturbed porous soil underneath the litter layer acts 
as a large storage area. Some of the water in this soil-storage area is 
available for use and transpiration by trees and other plant life. The 
excess water filters gradually through the soil and ultimately empties into 
streams which flow into the reservoir system. 

As the plantations reach 25 to 30 years of age, thinning or improve¬ 
ment cutting to maintain a thrifty stand of trees is practiced. However, 
before all mature trees are removed from a forest, stand, an adequate 
supply of natural reproduction has become established. Reinforcement 
plantings of seedlings obtained from a forest nursery are also used. Timber 
operations remove from the vatersbed such important wood products as pulpwood, 
firewood, poles for piling, and savlogs. 

A continuous development of the watershed management program is 
proceeding in a number of areas such as: (l) work toward a master plan for 
the sustained management of the watershed forests; (2) effective control of 
timber-harvest and timber-culture operations; (3) inspection and protection 
of property boundaries and watershed forests; (4) research and study of 
vegetative influences on the watershed; and (5) supplemental support of 
administrative and supervisory functions. 
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The role of forest cover in providing an adequate supply of high 
quality water is a major consideration of the Water Supply Board. There¬ 
fore, modern and efficient techniques are continually being applied in 
forestry and maintenance operation on the watershed. 

6. Intake and Treatment Works 


Water is conveyed from Scituate Reservoir by gravity and through 
aqueducts to the Water Treatment Works, where ferric sulphate is added 
as a coagulant. The chemicalT treated water is then subjected to an 
influent aerator to remove carbon dioxide which has corrosive properties 
and other gases which may produce disagreeable taste and odors. Influent 
aeration is also practiced to oxidize iron and manganese to be removed in 
the coagulation and sedimentation processes which follow. 

From the influent aerators, the water continues under gravity 
to a lPrge, circular mixer. The tangentially entering stream can be 
regulated to produce desired velocities to insure thorough mixing of the 
chemicals. Slaked quicklime is introduced to the water just ahead of the 
mixer. The lime aids in reducing the corrosive properties of the water 
by raising the pH of the water from an acid to an alkaline state. This 
treatment is necessary for the removal of iron and manganese which 
cannot be removed at a low pH, and affords a better degree of coagulation 
by increasing the specific gravity of the ferric hydroxide floe. 

After a detention period of two to three days in the coagulation 
and sedimentation basins, which have a combined capacity of 160.21 MG, 
the water is drawn through rapid sand filters, chlorinated and treated with 
sodium silico-fluoride, and then discharged into a clear well for 
distribution through the h.5 n ile aqueduct to the city. 

In case of breakdown or repairs to the purification works or the 
main aqueduct, there is reserve storage in three underground concrete 
reservoirs within the system. Neutaconkanut Reservoir in Johnston stores 
about 38.58 MG at its normal operating level; while Aqueduct Reservoir 
in Cranston holds ^0.03 MG at normal operating level. The third reserve 
storage reservoir, Longview, located in North Providence, contains 11.9*+ MG 
at normal elevation. The combined storage capacity of all three reservoirs 
is sufficient storage for two days' supply based on the average daily 
consumption and represents one day reserve at the rate of the maximum day. 

C. Method of Evaluation for Water Contamination 


1. Basic Assumptions 

The basic assumptions and method of evaluation are essentially 
those presented previously (,l)(67). Radionuclides from fallout are 
assumed to mix homogeneously in a water supply reservoir following deposition 
on the water surface. Therefore, the concentration of a radionuclide in 
water may be determined from a knowledge of the nuclide surface density, 
the reservoir surface area and volume of the body of water. 

A second assumption is that the parent elements of all long-lived 
radionuclides of interest have already decayed to a negligible amount 
at H + 1 hour. The nuclide concentration estimates are calculated from 
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fallout intensity contours corrected to H + 1 hour. This correction means 
that after local fallout stops, essentially 24 hours a^ter detonation, the 
existing fallout patterns are traced back to a common time basis at one 
hour after weapon detonation by use of a typical ionization rate decay 
curve. The following relationship may then be used to obtain activity 
concentration from nuclide concentrations in water: 


A = X N 

c 


( 11 ) 


where: 


2 . 


A = the activity concentration in water (curies per unit volume) 
X - the radioactive decay constant (time - ^) 

= nuclide concentration in water (atoms per unit volume) 

General Procedure 


The contamination in water supplies may be derived from any one 
of four sources or any combination of these: 

a. Reservoir Supplies: 

(1) Direct contamination 

(2) Contamination from fe* er streams 

(3) Contamination from watershed runoff 

(4) Contaminated ground water inflow 

b. Stream Supplies: 

(1) Direct contamination 

(2) Contamination from runoff 

(3) Contaminated groundwater inflow 

The evaluation is performed in two stages. First, the effect 
from a single weapon is analyzed, -"hen by the principle of superposition, 
the combined effect of several weapons is determined. 

The general procedure may be outlined as follows: 

a. The fallout pattern for each assumed weaponage is superimposed 
over the area map with ground zero coinciding with the target point, and 
the downwind axis parallel to the prevailing wind direction. 

b. Working from the fallout model, the soluble nuclide surface 
density is evaluated and integrated over the surface area cf interest and 
the result divided by the total volume of water to obtain the nuclide 
concentration in water. The activity concentration, A, is then estimated 
from Equation (ll)., 
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3. Runoff 

Contamination due to runoff is based on the assumption that all 
the nuclides which are distributed over the watershed enter the water supply 
system in one form or another. As part of this research program investigations 
are currently underway to determine to what extent this assumption is valid 
or if the number of nuclides deposited on a surface and entering a water 
supply system is significantly reduced from 100 per cent due to factors 
such as vegetative and soil uptake. 

Runoff volume is calculated from the relationship ft = C i A; 
where C is an average runoff coefficient determined from rainfall and 
runoff records; i is the estimated maximum rainfall intensity; and A is 
the watershed surface area. 

Hie combined effect of direct surface and runoff contamination 
is obtained by computing a weighted average of the radioactive concentrations 
calculated for each case. 

D. Absorbed Dose 


During a nuclear attack, public water supply systems are subject to 
severe damages and the contents exposed to fallout contamination. Soluble 
radionuclides mix with the feeder streams of a watershed and tend to 
increase the radioactivity in the water. When this water is consumed 
it may constitute a major source of the internal radiation hazard. 

A number of mathematical models for estimating the absorbed dose 
from assimilation of radionuclides in body organs of humans have been 
developed. Estimates reported here are based on the Miller-Erown Model 
of Biological Uptake and were determined according to criteria presented 
in a previous report entitled, "Evaluation of Fallout Contamination of 
Water Supplies" (l). 

E. Results 


The concentrations of six biologically important radionuclides in the 
Providence water supply system at H + 1 hour following a 30 MT nuclear attack 
at Springfield are presented in Table VII. 

Assuming a standard intake of one liter pe.. person per day, the absorbed 
dose for total body organs for different starting times after detonation 
t Q , and ingestion periods t, has been summarized in Table VIII. 

F. Discussion of Results 


A summary of the water contamination levels and decontamination 
requirements for the six biologically important radiosotopes investigated 
is presented in Table IX. As may be seen, the activity levels from 
surface contamination of 1-131 and Ba-l^+O exceed 10~? pc/mi, the tentative 
emergency standard according to Bale (68). Peacetime continuous occupational 
exposure MPC (69) values for 1-131 and Ba-1^0 in drinking water are 
6xlO"5 (ic/mi and 8xl0 _1+ pc/mi, respectively. 
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Although the contamination levels calculated for the Providence, R. I. water 
supply are higher than those calculated for other municipal water supplies (l) 
the results agree reasonably veil. It car. be expected that a more refined 
analysis of contamination to the Providence water supply (including the 
effects of stream flow between reservoirs, sedimentation, and tne factors 
affecting runoff) will reduce the giver, estimates to a range within those 
previously calculated for other municipalities, especially since ion exchange 
and plant uptake phenomena are expected to decrease the level of water 
contamination. 
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TABLE VII 


Contamination of Providence, 

Rhode Island, water Supply 


System 

Isotope 

atom/iiter (xlO* 2 ) 

pc/ml 

W'estccnnaug Reservoir 

Sr-89 

3.36 

1.37x10" 2 


Sr-90 

5.91 

i.26xl0~ 4 


Eu-106 

2.80 

1 .66xl0~ 3 


1-131 

6.15 

1.66X10' 1 


Cs-137 

4.28 

7.63XIO' 5 


Ba-140 

6.21 

1.05xl0 _1 

Pcnaganset Reservoir 

Sr-89 

I.69 

6.97xlO’ 3 


Sr-90 

3.38 

7.19xlO -5 


Ru-106 

1.62 

9.63x10 -1t 


I-??l 

3.55 

9-57xlO" 2 


Cs-137 

2.35 

4.19xl0" 5 


3a-l40 

3.57 

6.03x3 0" 2 

Barden Rese^oir 

Sr-89 

2.32 

9,43xlO" 3 


Sr-90 

4.49 

9.55x10" 5 


Ru-106 

2.32 

1.26xl0" 3 


1-131 

1.66 

1.26x10" 1 


Cs-137 

3.22 

5.74xlO -5 


Ba-l4o 

4.72 

7.97x10' 2 

Moswensicut Reservoir 

Sr-89 

0.546 

2.22x10" 3 


Sr-90 

1.04 

2.21al0" 5 


Ru-106 

0.488 

2.90xl0" 4 


1-131 

1.07 

2.88x30" 2 


Cs-137 

0.767 

1.37x10“ 5 


3a-l40 

1.09 

1.84x10" 2 
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TABLE VH (Cont’d) 


Contaminat 

ion of Providence, 

Rhode Island, Water Supply 

System 

Isotope 

l2 

atom/liter(xl0' ) 

uc/ml 

Regulating Reservoir 

Sr-89 

Sr-90 

Pu-106 

1-131 

Cs-137 

Be-l4o 

2.33 

4.46 

2.11 

4.63 

3.27 

4.70 

9.47x10”“ 

9.46x10" 5 
1.25xl0" 3 
1.25X10* 1 
5.83xlD~ 5 

7.94x10 ” 2 

Seituate Reservoir 

Sr-39 

0.683 

2.78xl0” 3 


Sr-90 

1.29 

2.74x10' 5 


Ru-I06 

0.6i4 

3.64x10“ 4 


1-131 

1.35 

3.64x10"^ 


Cs-137 

0.959 

1.71>cl0" 5 


Ba-140 

1.36 

2.30x10” 2 

Entire System 

Sr-89 

1.21 

4.91x10 -3 


Sr-90 

1.47 

3.13xlO" 5 


Ru-106 

0.698 

4.15x10” 4 


1-131 

1.53 

4.12x10" 2 


Cs-137 

1.08 

1.92xlO" 5 


Ba-l4o 

1.54 

2,60xl0" 2 

Entire System 

Sr-89 

18.1 

7.34x10" 2 

Including Runoff 

Sr-90 

34.9 

7.42x10" 4 


Ru-106 

16.6 

9.86x10" 3 


1-131 

36.4 

9.8lxlO" x 


Cs-137 

25.2 

4.49x10' 4 


Ba-l40 

36.6 

6.l8xl0 -1 
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.TABLE VIII 


Internal Hazard of To 

tal Body from Ingestion 

of Contaminated 



Providence Water Supply (rens) 



t 

f 

Dose due to direct 

Dose due to the effect 


O 

L 

0 

contamination 

of runoff 

Isotope 

(days) 

(days) 

of reservoir 

included 



30 

0.188 

2.80 


1 

91 

1.11 

l6.6 

Sr-89 

7 

30 

0.115 

1.72 


( 

91 

0.939 

14.0 


ik 

30 

0.0536 

0.80 


91 

0.761 

11.3 


T 

30 

0.00316 

0.0750 



91 

0.030 

0.712 

Sr-90 

7 

30 

O.OOI98 

0.0471 


91 

0.0260 

O.618 


l4 

30 

O.OOO986 

0.0234 


91 

0.0219 

0.520 


T 

30 

0.00244 

0.0581 



91 

0.00956 

0.227 

Ru-106 


30 

0.00173 

0.412 


l 

91 

0.00879 

0.209 


14 

30 

0.000977 

0.0232 


91 

0.00789 

0.188 



30 

1.10 

26.2 


X 

91 

1.51 

35.9 

1-131 

7 

30 

0.545 

13.0 


1 

91 

O.898 

21.4 


lb 

30 

0.203 

4.84 


91 

0.490 

11.6 


*1 

30 

0.00313 

0.0731 


X 

91 

0.0252 

0.587 

Cs-137 

7 

30 

0.00201 

0.0469 


( 

91 

0.0222 

0.519 


lb 

30 

O.OOO989 

0.0231 


91 

0.0191 

0.446 


*1 

30 

0.0524 

1.24 


X 

91 

O.O996 

2.37 

Ba-140 

7 

30 

0.0291 

0.692 



91 

0.0713 

1.69 


14 

30 

0.0124 

0.294 


91 

0.0480 

1.14 




I 





r 



* 
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Summary 

Isotope 
Sr-89 
3r-90 
Ru-106 

1-131 

Cs-137 

Ba-l40 


IABLE IX 


of Water Contamination Levels for the Providence, R.I., Water Supply 

(all values given are in pc/mi) 


Surface Contaga nr. tion 

4.91 x ID" 3 
3.13 x 10’ 5 
4.15 x 10~ k 
4.12 x 10" 2 

1.92 x 10" 5 
2.60 x 10" 2 


Contamination, 
including runoff 

7.34 x 10" 2 
-4 

7.42 x 10 
9.86 x 10' 3 
9.81 x ID" 1 
4.49 x ID" 4 
6.18 x 10" 1 


• -.Mmmcmaiun 
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